Pedestrian flow through multiple bottlenecks 
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We investigate the dynamics of the evacuation process with multiple bottlenecks using the floor 
field model. To deal with this problem, we first focus on a part of the system and report its 
microscopic behavior. The system is controlled by parameters of inflow and the competitiveness of 
the pedestrians, and large inflow leads to a congested situation. Through simulations, the metastable 
state induced by conflicts of pedestrians is observed. The metastability is related to the phase 
transition from free flow to congestion. The critical condition of the transition is theoretically 
derived. 

In addition, we give simulation results of situations with multiple bottlenecks. They imply that 
local improvement of pedestrian flow sometimes adversely affects the total evacuation time, and 
that the total optimization of the system is not straightforward. 
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I. INTRODUCTION 

Dynamical behavior of crowds has attracted many 
physicists over the last decades for its nontrivial charac- 
teristics [HII|- The motion of pedestrians can be regarded 
as a problem of a many-body system of "self-driven" par- 
ticles. In order to investigate the collective phenomena 
of the system, many microscopic models have been de- 
veloped: the social force model H, the floor field (FF) 

model [Ml 

, the lattice gas model [l3| . etc. In addition 
to the simulations, much effort has also been devoted to 
experimental studies [THHH]. In this research field, the 
evacuation of crowds has been vigorously studied since 
it is of great importance to design buildings properly for 
the case of emergency, in the context of risk management 
. One remarkable phenomenon observed during evac- 
uation is the clustering of pedestrians at bottlenecks such 
as exits (arching). When more than one pedestrian tries 
to move to the same place, the conflict occurs, decreasing 
the outflow of pedestrians. The effect of the conflict is a 
dominant factor of total evacuation time; however, only 
a few theoretical analyses have been performed so far 
In previous studies, the evacuation process from a 
single room has mainly been investigated. However, in an 
actual emergency, the pedestrian flow experiences many 
bottlenecks and merges together toward the exit. To the 
authors' knowledge, no systematic approach to evacua- 
tion from complex buildings has been proposed. This 
study provides a first step toward the understanding of 
the problem. First, as illustrated in Fig. [1] we abstract 
the most important factors, i.e., bottlenecks and their 
connectivity as a network. By considering this kind of 
a general network, one can apply our results to a broad 



range of practical problems. In this study, different from 
most other studies of networks, each node itself has a 
complex dynamics; we therefore begin by focusing on a 
single segment of the network. 
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FIG. 1. Schematic diagram of pedestrian evacuation from 
a building. A single segment enclosed with dotted lines is 
mainly focused on. 

In this paper, the overall argument is based on the 
FF model, which is one of the well-established cellular 
automata models for describing the pedestrian dynamics. 
The model is convenient not only for simplicity and ease 
of use, but for its extensibility |5l-fl2| . 

The effect of the conflict was first implemented in the 
FF model by introducing the friction parameter in Refs. 
[f| Q . In the present study, we use a more general ver- 
sion of the friction parameter, namely, the friction func- 
tion which was first proposed in Ref. Q. The motion 
of pedestrians involved in the conflict is canceled with 
a certain probability determined by the friction func- 
tion, which controls the strength of clogging and sticking 
among pedestrians. 
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In addition to exits, we set an entrance providing the 
system with pedestrians with a certain probability every 
time step. The stochastic entrance is convenient for us to 
control the inflow of pedestrians 0, [l(| , and by regarding 
it as inflow from an exit of the previous bottleneck, we 
can evaluate the effect of connectivity of the bottlenecks. 

The rest of this paper is organized as follows. Section 
im gives the definition of the model. In Sec. IIII1 sim- 
ulation results are shown. To explain the phenomenon 
analytically, we propose the second-order cluster approx- 
imation in Sec. IIV1 Finally, we summarize the argument 
in Sec. H 



II. MODEL 
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A. Floor field model 



We consider a two-dimensional lattice representing a 
room with an entrance and exit, consisting of N x N sites 
labeled = 1,2, •••,7V). Each site can contain 

only one pedestrian at most. Every time step pedestrians 
choose one destination site out of their five neighboring 
sites including the present site : (i + (i — 

1,J')> (iJ + !)> and (iJ - !) ( see Fi S- 0, according 
to two types of FFs. One of the FFs is the static FF 
(Sij) describing the shortest distance to the exit site, and 
the other is the dynamic FF (Dy ) expressing the total 
number of pedestrians who have visited the site. The 
dynamic FF has the dynamics of diffusion and decay, 
unlike the static FF p]. The transition probability 
for a jump to the neighboring site is determined by 
the following expression: 

Pij = Z^j exp {-k s Sij + k d Dij), (1) 

where k s and kd are non-negative sensitivity parameters, 
and Z stands for the normalization factor. £ returns 
for forbidden transitions such as to a wall, an obstacle, 
and neighboring occupied sites, and returns 1 for other 
transitions. In this paper, the static FF is given by the 
L 2 norm as 



Sij = \J [xij - x cx ) 2 + (yij - y ox ) 2 , (2) 

where (xij,yij) and (x ex ,y ex ) are the coordinates of the 
site and the exit site, respectively. On the other 

hand, we ignore the effect of the dynamic FF (kd = 0) 
since it does not greatly affect the arguments in this work 

m- 

B. Conflict resolution 

Due to the use of parallel update it happens that more 
than one pedestrian tries to choose the same site, which 
is called the conflict. The simplest solution of the con- 
flict is to choose one pedestrian randomly to move to the 



FIG. 2. Update rules. Each pedestrian can hop to its neigh- 
boring sites or stay at its present site in a time step. Pedestri- 
ans enter the area from the entrance with the probability a, 
and leave the area from the exit with the probability P — 1. 

site, and keep other pedestrians at their present sites. 
However, in actual situations, pedestrian flow is often 
clogged by more than one pedestrian moving at the same 
time. To model this effect, the friction parameter has 
been introduced 0, @ , and many significant results have 
been obtained so far. In a recent study Q, the friction 
function has been proposed to describe the effect more 
precisely. In the friction function, the number of pedestri- 
ans involved in the conflict, k, is reflected to its resolution 
probability. In this paper, we assume <p in the following 
form as in Ref . [1, Q : 

d>(<;,k) = i-(i-() k -k((i-c) k - 1 . (3) 

Here, £ G [0, 1] is the friction coefficient representing the 
strength of the clogging irrelevant to k. This is a mono- 
tonically increasing function of k and £. Note that this 
choice of 4> is one of the possible expressions. If one takes 
the friction function independent of k, it coincides with 
the friction parameter. 

Each conflict is resolved with probability 1 — (f>(£, k), 
and one of k pedestrians is randomly selected to move to 
the site, otherwise, the conflict remains. 



C. Entrance and exit 

In each time step, a pedestrian is provided to an en- 
trance site with the probability a 6 [0,1] if the site is 
empty, and removed from an exit site with probability 
(i. (See Fig. [51) In this paper, we assume that each 
room has enough space so that pedestrians at exit sites 
are smoothly accepted to the next room, namely, (3 — 1. 
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FIG. 3. (Color online.) Fundamental diagram of the system in the competitive parameter regime [(a)£ = 0.8] and the 
cooperative parameter regime [(b)£ = 0.0]. Each plot is obtained by setting the inflow probability a = O.f , 0.2, • • • , 0.6 and 
averaging the flux and density over 100 time steps. 




FIG. 4. (Color online.) Snapshots of the simulations. A 
green site (at the top) and a red site (at the bottom) indicate 
the entrance and exit sites, respectively. With the same inflow 
rate, one can observe the free-flow phase (left) and congestion 
phase (right) of pedestrians represented by blue circles. 

III. SIMULATIONS 

In the following, we set k s — 10 [l8j . The dimensions of 
the simulation area are 25 x 25, and one entrance site and 
one exit site are set at (13,25) and (13, 1), respectively. 
In this section we see some simulation results, varying 
the inflow probability a and the conflict coefficient 

A. Metastability of pedestrian flow 

In Fig. |31 the average pedestrian outflow through the 
exit, q, is plotted. Here, the density p is defined as the 
number of pedestrians in the room divided by the num- 
ber of sites in the area, 25 x 25. We performed simu- 
lations for 100000 time steps for each inflow probability 
a = 0.1, 0.2, ■■•,0.6 with the initial condition that no 
pedestrian is set in the system. The simulation results 
of the flux and density are averaged over every 100 time 
steps. The relation between the flux and density is often 
referred to as the fundamental diagram in the context of 
traffic flow, and in vehicular traffic the metastable state 
is observed in the fundamental diagram. The metastable 
state indicates an unstable state with high flux and den- 



sity before falling to the jammed state, and at the same 
density we can see multiple fluxes corresponding to the 
metastable state and jammed state, in a certain regime 
of density. Interestingly, the metastability is observed 
also in this problem. While in vehicular traffic flow, the 
metastability comes from the effect of inertia of vehicles; 
in this case, the conflict plays an important role for pedes- 
trian traffic. If the conflict occurs at the exit, outflow is 
decreased; therefore, even for the small density of pedes- 
trians in the system, the jammed flux can be observed 
by the concentration of pedestrians at the exit. On the 
other hand, if the pedestrians are dispersed and in the 
state of free flow, the system can keep large flux. This is 
the mechanism by which the metastability is induced in 
the pedestrian bottleneck . Consistent with these argu- 
ments, the metastability cannot be observed when ( = 
[Fig. Mb)}- 

Next, let us explain the fundamental behavior of the 
system. First, the free-flow phase (see Fig. 0} can be 
observed for every a and £. In this phase, large a directly 
leads to large flux, corresponding to the linear relation 
between q and p. On the other hand, when the conflict 
occurs at the exit, the rate of outflow shrinks. If this 
outflow is smaller than inflow (for large a and £), the 
density of pedestrians increases leading to congestion [the 
congestion phase shown in Fig. S (right)]. In contrast, 
if the inflow is not enough large, the conflict disappears, 
recovering the free flow. 

Here, if we take a as a controllable parameter for a 
given C, how can we optimize the inflow parameter? The 
answer is to keep the inflow lower than the critical rate, 
which prevents the conflict at the exit from growing. We 
discuss this issue in the following sections in detail. 



B. Critical phenomenon related to a and £ 

In this subsection, we analyze the phase transition 
from free flow to congestion in detail. Figure [S] shows 
the density and flux in the steady state and the transient 
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FIG. 5. (Color online.) Pedestrian flux and density in the steady state (st.) and transient state (tr.) vs a. Equations (5) and 
(12) derived in Sec. HVI are also shown as g/ theo. and q c theo., respectively. The curves are added to improve visibility in (b). 
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FIG. 6. (Color online.) Phase diagram obtained by simu- 
lations (st.). The red line indicates the theoretical critical 
condition (|13|) , 



state. Here, the transient state is denned as the state 
of the system with free-flow initial conditions (the room 
with no pedestrian) at finite time steps and introduced 
to include the probability of the system having the phase 
transition. Even in the parameter regime of the conges- 
tion the system might keep the metastable state for finite 
time steps. Hence the expected value of the flux (den- 
sity) in the transient state is larger (lower) than that in 
the steady state. These kinds of quantities are also im- 
portant because the system does not always reach the 
stationary state in actual evacuation of crowds. 

Here, we summarize the simulation conditions in Table 
|TJ In the simulations, each plot of the transition state 
(tr.) is obtained by averaging the flux or density over 
t m ax — 100000 time steps and 100 samples [l9j]. On the 
other hand, plots for the steady state (st.) are calculated 
by averaging over 1000000 time steps from t — 100000 to 



i max = 1100000. For the st. plots, we adopted the initial 
condition that pedestrians occupy all the available sites 
to ensure the occurrence of the congested situation in the 
corresponding parameter regime. 

In Fig. [5] (b), the lines of the density jump at each 
critical a, which correspond to the occurrence of the con- 
gestion. Corresponding to the fact, the flux of pedestri- 
ans [Fig. [5] (a)] agrees with two types of lines: In the 
free-flow phase, the flux is determined only by the in- 
flow probability; on the other hand, in the parameter 
regime of the congestion, the flux is determined by the 
outflow. These two pedestrian fluxes are evaluated the- 
oretically in the next section. Moreover, the simulation 
results of the steady states imply that the situation near 
the exit can be assumed to be independent of the in- 
flow probability in the congested situation. On the other 
hand, one can see the peak of the flux in the transient 
state. Since the system does not always cause congestion 
even in the congestion regime of parameters (metasta- 
bility), the expected flux is higher than the congestion 
flux for a around the critical value. If the inflow rate 
is large enough, the probability of the flux falling to the 
congested situation increases, and thus, the average flux 
decreases. 

Figure [S] shows the phase diagram of the system. For 
each C, the upper limit of the inflow rate to keep the free 
flow is depicted with the line derived in Sec. IIV CI From 
a practical standpoint, it corresponds to the criteria for 
avoiding the congested situation. 



TABLE I. Simulation conditions. 
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tr. 


averaging time steps 


4= 100000 - 1100000 


t = 


1 - 100000 


number of samples 


1 




100 


initial condition 


p = l 
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FIG. 7. (Color online.) Snapshots of the merging, (a) Center 
(ce) connection, (b) corner (co) connection. 



inflow is reduced by a strong bottleneck, the total outflow 
is improved conversely. 

Thus, we can conclude that pedestrian flow should be 
dispersed not only spatially but also temporally, and the 
strong bottlenecks might be used for the total optimiza- 
tion. Note that here we considered the situations that the 
effective inflow rate to the entrance hall is larger than the 
exit capacity, namely, they are in the congestion phase 
as a whole. If the exit has enough width, the local opti- 
mization directly improves the total evacuation time. 

IV. THEORETICAL ANALYSES 

Let us analyze the pedestrian outflow in each phase 
and its critical conditions focusing on a single bottleneck. 
In this section, we consider the limit k s —> oo, where 
pedestrians surely move to the most desirable site. 

A. Free-flow phase 



C. Merging 

Buildings have branches of passages, and the merg- 
ing of the pedestrian flow is an important and complex 
problem in architectural design. In this subsection, we 
present some simulation results to show a paradoxical 
effect of local improvement of the pedestrian flow. As 
more macroscopic systems than those in previous sub- 
sections, we consider the systems with three rooms and 
two types of connections: the center (ce) connection [Fig. 
[7(a)] and the corner (co) connection [Fig. [7(b)] . Here, 
two rooms are connected to a room with an exit site like 
an entrance hall, and we do not provide pedestrians to 
each room after setting initial conditions. In Ref. Q, 
it has been demonstrated that the exits at the corner 
has larger capacity of outflow than ones in the center be- 
cause they have only two neighboring sites, and conflicts 
are reduced. Therefore, at first sight, the corner connec- 
tion seems to be better in the sense of swift evacuation. 
In fact, however, it is quite opposite, as shown in Fig. 
ISJ The simulations have been performed with an initial 
condition that 50 pedestrians are randomly distributed 
in each of two rooms other than the entrance hall. As 
expected, the evacuation time from each room to the en- 
trance hall is improved by setting the connection at the 
corner, namely, t co < t ce ; however, the total evacuation 
time worsens as ( increases. Accordingly, the rate of the 
local evacuation time t co /t ce decreases while that of the 
total evacuation time T co /T ce increases. 

This phenomenon is explained as follows. The total 
evacuation time is determined only by the outflow from 
the exit, and the local pedestrian flow to the entrance hall 
does not influence the total outflow directly. Further- 
more, the large inflow to the entrance hall will increase 
the density of pedestrians near the exit, which leads to 
the decrease of the outflow. In other words, if the local 



In the free-flow phase, the pedestrian flux is evaluated 
by the inflow. First, we consider a balance equation at 
the entrance site: 

a(l - /9cn) = Pen- (4) 

Here, p en is the probability of finding a pedestrian at the 
entrance site. Since we assume k s — > oo, a pedestrian at 
the entrance site surely leaves the site (with probability 
1) in the next time step. By solving the equation, we can 
evaluate the flux qf\ 



This expression is shown in Fig. [5] (a) and well agrees 
with the simulation results. 



B. Congestion phase 

In the congested situation, the area near the exit is al- 
most fully occupied. This fact enables us to estimate the 
outflow in the steady state. In the previous studies, the 
pedestrian density at the neighboring sites of the exit is 
approximately assumed to be 1 d-Q. However, when the 
effect of the conflict is strong, this assumption does not 
give a good estimate. Especially in this paper, since the 
effect of the conflict depends on the number of pedestri- 
ans involved, we have to take the configurations of pedes- 
trians at the exit into consideration. Hence, we adopt 
a second-order approximation here. In the approxima- 
tion, we assume the probability of finding a pedestrian 
at the site where the exit is accessible with two jumps 
(the second neighboring site) to be 1 (see Fig. [S]). Then, 
the states of the neighboring sites are characterized by 
four occupation numbers A, B, C, and E, which take 
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FIG. 8. (Color online.) (a) Average evacuation time and (b) rates of evacuation time. T and t are the average total and local 
evacuation time, and the subscripts "co" and "ce" stand for the conditions of the corner connection and the center connection, 
respectively. Each plot is calculated by averaging 1000 samples. 



(empty) or 1 (occupied). By considering transition prob- 
abilities among these states, we can obtain the probabil- 
ity distribution P^bc °f finding each state in the steady 
state. Here, the superscript corresponds to the occupa- 
tion number of the exit, and the subscripts indicate one 
of its neighboring sites as shown in Fig. [S] To reduce the 
dimension of the transition matrix, we use the following 
facts. First, we can easily find Pq Q0 = and Pin = 0. 
Since no configuration of pedestrians can result in these 
states in the next time step, they are not realized in the 
stationary state. Furthermore, by symmetry of the sys- 



tem, the equations 
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are satisfied. With the normalization condition, the tran- 
sition matrix is summarized as follows: 
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Note that abbreviated notations 02 = 0(C, 2), 03 = 0(£, 3), 02 = 1 — 02, and 03 = 1 — 03 are used. Then, the pedestrian 
flux is given by 

I 

hand, when the effect of the conflict is not strong, the 
approximation gives good evaluation. 



1c — Pooo + 2-P 100 + P 010 + 2P 110 + P 101 . 



(12) 



This expression is illustrated with simulation results in 
Fig. [5] (a). When the effect of the conflict is dominant, 
the error becomes large. In this parameter region, pedes- 
trians are not provided smoothly to the second neighbor- 
ing sites due to the strong friction, and thus, the assump- 
tion of the approximation is not satisfied. On the other 
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FIG. 9. Cluster approximation. The probability of finding 
a pedestrian on the second neighboring sites of the exit is 
assumed to be 1. 



C. Critical conditions 

In the previous subsections, we evaluated pedestrian 
flux in the free-flow situation qf and congested situation 
q c . When q c becomes smaller than qf, the system cannot 
maintain the free-flow situation; therefore we can obtain 
the critical condition regarding the inflow probability for 
a given £ from the condition, q c = q/' 

a CI = -^. (13) 

1 - q c 

This critical condition is compared with simulation re- 
sults in Fig. |51 and well describes the phase transition. 

Using this expression, we can control the inflow, pre- 
venting the occurrence of the congestion. For example, 



by securing additional routes and keeping the total in- 
flow to one bottleneck lower than the outflow in the 
congested situation, the clustering at the exits can be 
avoided, which may shorten the total evacuation time. 



V. CONCLUSIONS 

In this paper, we have taken a first step to treat the 
problems of multiple bottlenecks in evacuation process. 
To consider the problems, we focused on a part of the 
bottlenecks, using a stochastic entrance. Through sim- 
ulations, the metastable state of pedestrian flow arising 
from the effect of pedestrian conflicts is demonstrated. 
Supported by approximate analyses we have derived the 
expressions for the pedestrian flux in the free-flow phase 
and congestion phase, and a critical condition of the in- 
flow to prevent the congestion. Furthermore, interesting 
phenomena related to the merging of pedestrian flow have 
been reported. The local improvement of the pedestrian 
flow sometimes causes more serious congestion. We be- 
lieve we could give some hints for designing evacuation 
routes from a theoretical point of view. To validate these 
characteristics, some experimental studies would also be 
necessary. 

In this paper, only simple situations have been focused 
on to concentrate on the essence of the problem. How- 
ever, other factors such as the width of exits, multiple 
entrances and/or exits, route choice of pedestrians, ob- 
stacles, etc., should also be investigated in combination 
with our study in future works. 



APPENDIX: EXPLICIT EXPRESSION OF q c FOR A SIMPLE CASE 

Although we can obtain the explicit expression of q c by solving Eqs. (|10p and the expression is very complicated 
in general. If we use the friction parameter \x 0-0| independent of the pedestrian number, namely, <p2 — 4>z — /i, the 
expression will be relatively simple. Here we show the result in the second-order approximation for this case, <?2(/Uj 
and compare it with the previous result derived by first-order approximation [6J. From simultaneous Eqs. (|10[) and 
(fTTj). we can derive 

48 + 72/z - 132/i 2 - 28// + 140// - 236// + 131// + 49// - 91// + 125// - 126/i 10 + 57/i 11 - 9// 2 

q2 ^ P ' ~ 96 + 192/i - 144// - 68/i 3 + 240// - 404/i 5 + 78// + 129// - 166// + 185// - 117/i 10 + 48/i 11 - 9// 2 ' 

(14) 

On the other hand, Ref. [|| presented the expression derived by the first-order approximation, qi (/i) = ■ The 
difference between these equations yields 

1i M - <7i(m) 

//(32 + 16/t - 84// + 64// - 10// - 75/i 5 + 75/i 6 - 18//) 

~ (2 - /t)(96 + 192/i - 144/i 2 - 68// + 240// - 404// + 78// + 129// - 166// + 185// - 117//° + 48// 1 - 9// 2 ) 

= \^-\^ + --- ^ « l ) ( 15 ) 
• ' ' (16) 
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Hence (feO-t) = 33^ + \^ + 0{^ & ) is asymptotic to qi(fj) in the fi — > limit, and we can conclude that the second-order 
approximation presented in this paper improves the prediction, especially when the effect of the conflicts is strong. 
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